Characterization of the stretched exponential trap-time distributions in 
one-dimensional coupled map lattices 
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Stretched exponential distributions and relaxation responses are encountered in a wide range of 
physical systems such as glasses, polymers and spin glasses. As found recently, this type of behavior 
occurs also for the distribution function of certain trap time in a number of coupled dynamical 
systems. We analyze a one-dimensional mathematical model of coupled chaotic oscillators which 
reproduces an experimental set-up of coupled diode-resonators and identify the necessary ingredients 
for stretched exponential distributions. 
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I. INTRODUCTION 

The decay of certain quantities characteristic of man: 
complex systems such as glasses spin glasses 
quasi-crystals j|] , trapping models |[ _^ , coupled 



11 and 



non-linear systems [Qi S |9|' turbulence |1£ 

others [|l^ ^ Q is often described by a stretched- 
exponential - or Kohlrausch - functional form: 



(bit) ^ e^p[^it/Tf] , <>0 



(1) 



with < /3 < 1 Q. Although Eq. (Q) provides a good fit 
to a wide range of experimental and numerical results, in 
many cases these can also be fitted by power laws with 
comparable accuracy; while most experimental set-ups 
can span many decades in time, few have achieved more 
than 2 or 3 decades in (j){t). The Kohlrausch form can 
also be reproduced explicitely by a few theoretical models 
(see, e.g., Refs. |[ ^, ||, |l5| |6|) on the basis of various 
assumptions, but it is still unclear whether there can be 
a unifying solid theoretical justification for it. 

Recently, Hunt, Gade and Mousseau 0] found that 
stretched exponentials could fit experimental distribu- 
tions of trap time — the time a system spends in an 
uninterrupted state with temporal period two, in a one- 
dimensional network of coupled diode resonators — over 
more than 6 decades in the distribution. Numerical mod- 
els, based on the set-up, could expand this fitting over 
about 10 orders of magnitude, strongly suggesting that 
dynamics with underlying stretched exponential distribu- 
tions could be universal in coupled chaotic systems and 
ruling out a power-law or any other standard fit. 

In this paper, we revisit one such numerical model and 
provide a detailed characterization of the dynamics of 
this network as a function of system size and parame- 
ters, providing elements of explanation regarding the ori- 
gin of stretched exponential distributions in this system. 
In particular, we show that (1) size effects are important 
only for relatively small systems; (2) the natural invari- 
ant density, p{x), generated by typical orbits has well- 
defined structure that is self-organized; (3) the structure 



of p{x) alone cannot lead to stretched-exponential distri- 
butions — the dynamical spatial organization is essential 
for stabilizing the periodic orbits. 

This paper is organized as follows: In the next section, 
we introduce the mathematical model: a one-dimensional 



coupled map lattice. Section III discusses the background 
for the trap-time distributions. In section IV, the results 



from the simulation are presented, to be discussed in sec- 
tion |^. We summarize our conclusions in Sect. 



VI 



II. MODEL 

Our model is a one-dimensional chain of N diffusively 
coupled nonlinear deterministic maps, f{x), with a cou- 
pling constant a and periodic boundary conditions. The 
interaction is totalistic and involves only the nearest 
neighbors. The time evolution of this system is discrete 
and is described by an iterative equation: 

d 

Xn{t+1) - {l-a)f[Xn{t)] + - {/K_i(t)] + /K+i(t)]} . 

(2) 

which, given initial conditions x„(0) for each site n = 
1,2,...,A^, generates a time series {xn{t)}, for integer 
values of the time index t = 0,1, . . .. This model offers 
a much simplified version of the ID array of diode res- 
onators studied in Refs. |l^, ^ |l|, Although the 
individual diode resonators are best described by an in- 
ertial equation [ pT| , it was shown recently that Eq. 
captures the dynamics in the regime of interest here [[? 

The diffusive nature of the coupling in Eq. (|^) becomes 
more apparent if this equation is rewritten in the follow- 
ing form 

Xnit+1) = f[Xn{t)] + 

I {f[Xn-l{t)] - 2/K(i)] + fiXn+lim 



f[xn{t)] + -D+D^f[x,M 



(3) 



where the central difference operator Z3+Z?_/(x„) is the 
discrete Laplacian of f{xn) on a one-dimensional grid 



2 



with unit spacing. In the hmit N oo the discrete 
Laplacian can be substituted by its continuous coun- 
terpart and Eq. becomes a nonhnear analogue of a 
difference-differential diffusion equation, see e.g., Ref. |2^. 
For a large N, 1/r can be thought of as a "viscosity" and 
a as a "diffusivity" with their usual relation to the "tem- 
perature." 

We use the logistic map f{x) 



(a) 
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rx{l — x), Ref. 
order to describe the dynamics of the basic chaotic el- 
ements. Although this differs from the form studied in 
Ref. 1^, = 1 — ay^, the results are unaffected; the two 
maps are conjugate, related by a simple algebraic trans- 
formation: X = y/ay/ a/t -I- 1/2, g{y) = f[x{y)] - r/4 -|- 1. 
We choose the logistic map because it is more studied in 
the literature. In the simulations presented here, we fix 
the value of the coupling constant to a = 0.25, and vary 
N and r, which we refer to as the nonlinearity parameter 
in the following. 

Initial values of a;„(0) are taken from a random distri- 
bution in the [0, 1[ interval. Runs are then iterated for a 
few hundred thousands steps in order to avoid any tran- 
sient effect before starting the accumulation of data. All 
simulations presented here are done on one-dimensional 
arrays with variable length and periodic-boundary con- 
ditions. Statistics are generally accumulated over 10 mil- 
lion to 10 billion time steps. 

Following the experiment 0, the analysis of the dy- 
namics is done using coarse-grained variables defined by 



cr„(t) = sign[Xnit) - Xthr], t = 0, 2, 4, . 



(4) 



where the quantity x„(i) is defined in Eq. (g) and Xthr 
is a certain threshold value. The results presented here 
are not very sensitive to the value of the threshold. For 
simplicity, we select the value of the unstable fixed point, 
X* — 1 — 1/r, for the single map. This coarse-graining 
reduces the problem from continuous to two-state. 

The basic dynamics of the coupled oscillator being pe- 
riod two, the analysis is also done only over even (or odd) 
time steps. 



III. BACKGROUND: THE DISTRIBUTION OF 
TRAP TIME 
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FIG. 1: Coarse-grained time evolution of the coupled logis- 
tic map. The coarse-grained variable a„{t) in a stroboscopic 
representation : for every other iteration t a dot at (n, t) 
corresponds to a„{t) = -1-1 and the blank space signifies 
(Tnit) = -1. iV = 1000 and (a) r = 3.83, (b) r = 3.i 



of the microscopic density distribution of particles in a 
material psj . 

The relation between the auto-correlation function 
C{t) — {a{t')a{t' + t))t' of a renewal process a{t) — 
defined analogously to Eq. — and the distribution of 
trap time of this process was recently studied for a range 
of distributions ||2^. Using this framework, it is possible 
to show that C(t) corresponding to the stretched expo- 
nential distribution of occupation time is also stretched 
exponential at long time albeit with a different stretching 
exponent (3. 

These results are only valid when the traps are uncor- 
related in time which is the case for the systems studied 
here. The details of this work will be presented else- 
where [ p6| . 

This establishes a direct relation between the distribu- 
tion of trap time, discussed in this paper, and the more 
standard auto-correlation function, measured in a num- 
ber of experiments on glasses and other complex systems. 



We are interested in the statistical distribution of trap 
time in a coarse-grained state space of a one-dimensional 
chain of coupled nonlinear maps, as shown in Fig. ^ This 
quantity is formally equivalent to the distribution of time 
intervals between zero crossings of renewal processes such 
as random walks |24 and has the advantage that it can 
be measured experimentally and numerically to a high 
degree of accuracy for this system Q . 

Generally, however, experimentally measurable relax- 
ation responses are mathematically described by auto- 
correlation functions of certain dynamical variables. For 
example, the inverse Fourier transform of the dynam- 
ical structure factor, the intermediate scattering func- 
tion, can be calculated as time auto-correlation function 



IV. RESULTS 

Before discussing the possible origins for the stretched 
exponential distributions in our model, it is necessary to 
offer a characterization of its dynamics as a function of 
the parameters of the model. 

First, we consider the size effects on the dynamics 
of the network, studying experimental set-up varying in 
length between 15 and 256 oscillators, and simulated ar- 
rays of up to tens of thousands of sites. 

Then we investigate the dynamics of the system as a 
function of the nonlinearity parameter, r, which brings 
the system through a series of dynamical changes from 
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FIG. 2: Spatial correlation function C{n) — 
{(yn'{t)<^n+n'{'t))n' ,t- Dots corrcspond to the numerical 
data. Solid line interpolates numerical points and is shown 
as a guide to the eye. Dashed lines show the exponential 
decay of oscillations: ±0.55 exp(— n/8.25). Inset also shows 
C(n), but calculated for a different value of r. The am- 
plitude of oscillations (dashed lines) in this case decays as 
±0.5exp(-n/1.8). 

stable periodic orbits to full chaos. In particular, it is im- 
portant to assess the parts of the parameter space where 
the stretched exponential distributions can be observed. 

Once the basic phase diagram is established, we discuss 
the building up of the dynamics on a single site embedded 
in the network as well as the spatial structure associated 
with the stretched exponential distribution. 



A. Dependence on A'^ 

As is seen in Fig. |l|, the traps are directly associated 
with a periodic spatial organization, which is controlled 
by the coupling a. 

As a is increased, the spatial organization goes through 
two rapid transitions, showing qualitatively different fea- 
tures: For a < 0.1, the lattice tends to follow a period- 
two spatial organization. For intermediate values of 
0.1 < a < 0.19, the lattice immediately freezes into a 
period-two state in both space and time, for all values of 
the driving parameter, r. Above this threshold, a > 0.19, 
the dynamics becomes stochastic again while dominant 
spatial period goes to 4 and even longer for large a. In 
each of these phases, the variation of a affects only min- 
imally the spatial structure of the phase. 

In spite of the evident organization seen in Fig. g, the 
spatial correlation is short-range. Figure ^ shows that 
spatial correlations vanish exponentially fast with a typ- 
ical length scale between about 2 and 8, i.e., the only 
static spatial correlation appearing in the system is di- 
rectly associated with the short-range organization. 
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FIG. 3: Bifurcation diagrams for (a) an isolated map, (b) 
iV = 8, (c) iV = 16, (d) N = 32, (e) = 1024, and (f) 
= 10000. Dashed line in all panels shows the unstable 
orbit X, — 1 ~ 1/r. 



This sort-range correlation implies that size effects 
should be very limited. Fig. |^ shows the lattice-size de- 
pendence on the bifurcation diagram, for a single site on 
the lattice. 

For small lattices, 1 < iV < 32, the dynamics depends 
on whether N is even or odd. For even N, arrays syn- 
chronize rapidly while the odd sizes continue to display 
chaotic trajectories. For example, the bifurcation dia- 
gram changes qualitatively as one goes from a single iso- 
lated site to a chain of 8 or 16 oscillators: the chaotic 
phase disappears totally and the system remains peri- 
odic in the coarse-grained state space until r = 4.0 [see 
Fig. ^(b) and (c)]. While isolated oscillators produce 
an exponential trap-time distribution, this distribution 
tends to the stretched exponential form at short times 
for lattices with odd N as small as 15. Size effects are 
still present for these lattices, however, and the long-time 
distribution diverges from the stretched exponential [see 
Fig. ^ (a)]. Interestingly, the sign of the deviation for the 
stretched exponential oscillates as the array is increased 
in size two by two, indicating a certain spatial frustration 
in these small systems. 

For even N that is not a multiple of four, the mismatch 
with the periodicity results in the presence of inclusions 
of stable defects, with four sites at a row in the same 
band. These defects can also be present and stable in 
lattices where is a multiple of four. 

As the number of elements reaches 32, orbits corre- 
sponding to larger values of r^in < < 4 become chaotic 
and the bifurcation diagram approaches that of an in- 
finite lattice, rmin continues to decrease with growing 
number of sites, and is well converged for a lattice of a 
few hundred oscillators. As shown in Fig. ^ (b) the 
same trend is seen experimentally, although the inherent 
disorders helps to decrease the finite size effects. 

Even though the bifurcation diagrams for an isolated 
map and the infinite lattice appear similar, their respec- 
tive orbits are qualitatively different for most values of 
r. In particular, most of the chaotic region of the latter 
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FIG. 4: (a) Distribution of trap time for r = 3.8888. N = 21 
(squares: set of data points lying above the solid line) and 
A'^ = 19 (open circles: the set of data points lying below 
the solid line). The solid line corresponds to the stretched 
exponential fit Eq. ^ with /3 = 0.50 ±0.05 and r = 3.0 ±0.5. 
In inset, we show the data for A'^ = 17 (dash-dotted line), A'' = 
21 (dashed line) and the fit (solid line), (b) Distribution of 
trap time for the experimental set-up (squares: A'^ = 21, dots: 
A'' — 20, open circles: A'" = 19). The stretched exponential 
fit is with 13 = 0.75 ± 0.05 and r = 5.0 ± 0.5. Vertical lines 
in both panels demarcate the interval over which the fitting 
procedure was performed. 
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FIG. 5: Larger dots: bifurcation diagram of the map x{t+l) = 
0.75 f[x{t)] + 0.1625. Smaller dots: the same as in Fig. | (c). 
Although p{a) [dashed-dotted line in Fig. ^ (c)] corresponding 
to the bifurcation diagram depicted by the smaller dots is not 
a 5-function, it is narrow enough so that it creates a period- 
two attractor in the coarse grained phase space of the couple 
map chain. 
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shows a stretched exponential trap-time distribution. 



B. Dependence on r 



FIG. 6: Trap-time distributions fitted with stretched expo- 
nentials. Dots correspond to A'' = 1000, r — 3.83 and 
l3 = 0.33 ± 0.05, T = 11.5 ± 0.5. Open circles correspond 
to N = 64, r = 3.8888 and P = 0.50 ± 0.05, r = 2.9 ± 0.5. 
The solid lines show the best fit to the stretched exponential 
function ^(t). 



Fig. H also indicates the effect of r on the dynamics 
of a single site in a chain. The coupling stabihzes the 
orbits, reducing significantly the size of the chaotic re- 
gion. Its effect is to shift the bifurcation diagram to the 
right, moving the dynamics towards periodic orbits as is 
shown i n Fig . || (Fig. |^ will be discussed in more detail 
in Sect. |V^). 

The stretched exponential behavior appears, for a large 
enough lattice, around tq = 3.83. At this threshold value. 



the stretching exponent (3 w 0.33, obtained from fitting 
the long-time part of the trap distribution (Fig. |^) ac- 
cording to Eq. (|l|). 0.33 is the lowest value of (3 we could 
observe for the explored regions of the parameter space of 
the model. Below this threshold, chaotic windows are in- 
terspersed with periodic ones. Trajectories in these nar- 
row windows tend to fall erratically onto neighboring pe- 
riodic orbits, generating multiple-step trap distributions 
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with no clear overall time behavior. 

As seen in Fig. |[ the value of (3 increases with r. For 
r = 3.8888, the trap-time distribution is well fitted by a 
stretched exponential with (3 « 0.5, while at r = 4, the 
full chaos limit for an isolated logistic map, a fit to p{t) 
gives P = 0.70 ± 0.05 and t = 3.0 ± 0.5. 

Because of the short spatial correlation and the excel- 
lent quality of the simulation data, the exact value of /3 is 
well defined: increasing the lattice size from 64 to 1000 
sites leaves the trap-distribution essentially unchanged. 
Moreover, the same values of f3 are obtained, within the 
error bars, by changing the length of the interval over 
which the fitting procedure is performed. 

These numbers are also in reasonable agreement with 
experimental results which show P varying from 0.10 ± 
0.05 to 0.95 ± 0.05. This wider range for the experiment 
is probably caused by the presence of site disorder. 

This confirms that the best functional form for the 
simulation data is given by the Kohlrausch function. 

C. Single site dynamics. 

The iterative rule Eq. (|2|) can be viewed as an equa- 
tion describing the dynamics of a single element in the 
presence of an external additive perturbation a{t): 

x{t + l) = {l-a)f[x{t)]+aa{t) (5) 

where a{t) = {/[a;„_i(t)] + f[xn+i{t)]}/2. For simplic- 
ity, since all sites are statistically identical, we drop the 
subscript n in Eq. (||). 

This form allows us to concentrate on the impact, at 
the single-site level, of the rest of the network and to try 
to identify the essential elements for a stretched expo- 
nential dynamics. For this purpose, it is convenient to 
use the natural invariant density p(x) generated by typi- 
cal orbits {x{t)}, t = 1, 2, ... of the map in Eq. (|j). The 
notion of the natural invariant density is widely used in 
the studies of chaos |2^. The function p{x) is defined so 
that for any interval [x, x + dx] £ [0, 1] the fraction of 
the time typical orbits spend in this interval is p{x)dx. 
In the same way one can define the density p{a) for the 
perturbation a{t) in Eq. (^. 

Fig. shows the behavior of p{a) and p(x) for a few 
sets of the parameters a, N, and r. In the case of an 
isolated map (a = 0) and in the chaotic regime, a{t) is 
the sum of two identically distributed independent ran- 
dom variables — p{a) is shown by the dash-dotted line 
in Fig. (b). According to the central limit theorem, if 
this sum contained much more than two term, p{a) would 
very well agree with the normal (Gaussian) distribution. 
However, even with only two terms, p{a) exhibits a max- 
imum at a = 1/2 [the dash-dotted line in Fig. |^ (a)]. 

The presence of one or few peaks in p{a) makes the 
orbits {x(t)} intermittent for some values of the nonlin- 
earity parameter r. This intermittency means that the 
(coarse-grained) orbit {a(t)}, see Eq. (Q), stays periodic 
during some time, becomes chaotic for a certain while, 



and goes periodic again. In order to demonstrate how a 
periodic attractor develops in a coupled system at a value 
of r for which the isolated map produces a chaotic orbit, 
we can imagine a density p(a) = S{a — ao), where S{a) is 
the delta function and ao G [0, 1] is a constant. In this 
case, Eq. (||) becomes x{t + 1) = (1 — a)f[x{t)] + aao — 
cif[x(t)] +C2, which is the iterative equation for the iso- 
lated map scaled by a constant ci and shifted by another 
constant C2. This manipulation with the isolated map 
results in a bifurcation diagram translated with respect 
to the original one [Fig. I (a)] while preserving all its 
qualitative features, see Fig. H. In particular, the period- 
two attractors of the modified map overlap with chaotic 
orbits of the original map. 

As discussed above, in the interval 0.1 < a < 0.19, 
the lattice falls into period-two spatial and temporal 
patterns, i.e., each two neighboring sites oscillate out 
of phase and each site alternates between two values 
x^^'' < x^'^\ as indicated in Fig. |7| (a) and (b). 

For larger values of a, the natural invariant density 
becomes more complex. When the lattice size is a small 
even number, the spatial pattern is periodic with the pe- 
riod which is equal to four (two sites "up" followed by two 
sites "down" ) and each site again undergoes a period- two 
dynamics [Fig. ^ (c) and (d)] which is manifested by the 
peaks in p{x). For N = 16, there are narrow peaks in 
p(a) and p{x). These peak broaden and merge rapidly, 
however, for large N and the system displays strong sign 
of chaoticity. 

This chaoticity can be controlled by r and we see that 
the orbits spend more and more time in period-two at- 
tractors relative to the chaotic ones as the nonlinearity 
parameter is decreased, relating directly to a change in 
the value of the stretching exponent /3 [Fig. |^ (e) and (f )] . 

Although providing a clearer picture of the dynam- 
ics, knowledge of the invariant density is not sufficient to 
reproduce the stretched exponential dynamics. For ex- 
ample, choosing randomly the values of a{t) in Eq. (||) 
according to a prescribed p{a) generates only an exponen- 
tial trap-time distribution. Variations on this theme, in- 
cluding the introduction of a two-step distribution, which 
favors staying in the period- two phase once in it, also 
fail to give a stretched exponential distribution of trap 
times. Spatial correlations, even though short range, are 
essential to induce a stable and self-organized stretched 
exponential distribution. 



D. Stability of spatial period-4 structure. 

Fig. |l| plots the time sequence for cr„(i) on a small 
region of a 1000-site chain for r = 3.83 {(3 ~ 0.33) and 
r = 3.8888 {(3 ~ 0.5). These show a remarkable range 
in the trap size: the larger traps of the first one can be 
up to 250 times longer than those dX r — 3.8888, and 
corresponding to a time-scale 30,000 times larger than 
the basic time-step. This time scale is similar to that 
observed experimentally [Q . In spite of these very stable 
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FIG. 7: Densities p{a) and of a(t) — {f[xn-i{t)] + /[a;„+i(i)]}/2 [panels (a),(c),(e)] and x(t) [panels (b), (d), (f)] from 
Eq. (j^ normalized by the maximal value. Note that each row of panels corresponds to the same sets of parameters. 



traps, for long enough time, all sites display an identical 
dynamics and no region of the lattice remains frozen. 

Looking at the same figure, we also see that, for a 
given r, the trap width appears to be uncorrelated with 
its length: long traps remain narrow and are composed of 
ordered domains with several occurrences of period 4, 5 
or 6. While the period-4 domains can remain totally sta- 
ble for a long time, the middle site of the 3-up or 3-down 
segment in period 5 and period 6 domains show some 
instability. They must therefore be considered more as 
a defect in the stable phase than as an additional spa- 
tial structure. Even so, these larger basic domains, that 
occur mostly for the lower values of r, are responsible 
for the apparent spatial periodicity of 5 in the spatial 
correlation function at r = 3.83 (see Fig. 

In view of this discussion the space-time diagram can 
be separated in two distinct phases: the ordered spatial 
period 4 and the chaotic phase. The former can only be 
destroyed at the boundary or by spatial defects (domains 
of period 5 or 6). 

We verify this observation by two simulations. First, 
the entire chain is initially set into period-4 state with the 
period formed by two neighboring sites with a;„ (0) — x^^\ 
n = 1,2 and the next two sites with x„(0) = a;^^\ 
n = 3, 4. The values x'^^^ and x^^^ are taken to be equal 
to those corresponding to the two peaks in p{x) in Fig. |^ 
(f) for r = 3.83 and N = 1000. Simulations starting 
from this initial conditions remain frozen in this spatial 



period-4 and temporal period-2 state indefinitely. Al- 
though metastable, this state requires the presence of 
defects to be destroyed. 

This can be checked in a second test. We first iterate 
a 1000-site lattice at r = 3.83, for 10^ time-steps in or- 
der to eliminate any transient effects. At that point, at 
once, half of the chain sites are set into the same spatial 
period-4 state as above. As shown in Fig. ||, this phase 
then disappears gradually, over about 30,000 time steps, 
invaded from the edges by the chaotic phase. 



V. DISCUSSION 

The previous sections have shown that the dynamics 
of the network can be understood in terms of a competi- 
tion between two dynamical regimes: a stable period-two 
orbit and a fully chaotic state. This can be deduced from 
Figs. I and 0: when the neighboring sites are in opposite 
bands, their contribution shifts the bifurcation diagram 
into a period-two regime. A spatial period 4 is therefore 
the basic stable motif with this set of parameters. 

As mentioned above, we cannot reconstruct a stretched 
exponential dynamics simply using the natural invari- 
ant density inserted into Eq. (||). We have found that 
the only way to obtain a stretched exponential trap-time 
distribution in this situation is to impose a stretched ex- 
ponential trap-time distribution on a[t) for the values of 
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heterogeneities increases and the stretching exponent /3 
decreases under coohng in glass-forming systems. 

As we showed also, once in a perfect spatial period-4 
regime, the system will never become chaotic. The fi- 
nite life-time arises from defects in the period-4 phase 
or from the chaotic boundaries. This behavior is similar 
to what is seen in the spin-system models with effec- 
tive constraints on the dynamics and the state-space 
partition model Interestingly, both models show a 
stretched exponential decay of certain statistical quanti- 
ties. The distinctive feature of the present model. Sect. || 
is that the effective dynamical constraints arise in the 
course of the evolution described by the dynamical equa- 
tion (Q), unlike the above cited models where the fixed 
constrains are imposed on the neighboring sites and the 
dynamics is due to a Monte Carlo procedure. 



FIG. 8: Space-time diagram for iV = 1000 and r = 3.83. 
Dots correspond to trap boundaries — points in time when 
a site interrupts its temporal period-2 dynamics. Every fifth 
site and every 40th iteration are shown. The heavy horizontal 
line shows the region which was set into spatial period-4 state 
at a certain moment in time. Note that thus induced spatial 
period-4 domain does not have any bursts in the interior, and 
is affected by the chaotic phase only at the boundaries. 



a within the peak region in p(a) (see Fig. |^). More in- 
terestingly, connecting a single site, with a unidirectional 
coupling, to two sites selected at random on a lattice, is 
sufficient to induce a stretched exponential trap distribu- 
tion on this single site, albeit with a larger /? than that 
for the elements which belong to the lattice. 

Thus once an external perturbation, following a 
stretched exponential dynamics, is imposed on a chaotic 
element, the latter will immediately adopt a similar dy- 
namics. However, stretched exponential distributions 
cannot be observed in a self-organized process without 
spatial organization. For example, connecting the nodes 
of a balanced binary tree unidirectionally towards its 
root results in the densities p{x) and p{a) similar to 
those depicted in Fig. |^, (e) and (f) already for two lev- 
els in the tree. Nevertheless, the trap-time distribution 
measured at the root node is a pure exponential. This 
suggests that the stretched exponential distributions re- 
quires some spatial organization allowing time-limited 
spatial period-4 structures to occur and to be stabilized. 

The nature of this spatial organization is somewhat 
paradoxical. Correlations are very short range and the 
width of these domain is only weakly varying with the 
stretching parameter. While traps at r = 3.83 are up to 
250 times longer than those at r = 3.8888, their width 
is only 4 to 5 times larger, as displayed in Fig. |l|. This 
decoupling between spatial and temporal correlation is 
reminiscent of dynamical heterogeneities observed in in 



glass- forming systems [|16|, g7|, |28|. Associating a "tem- 
perature" with the control parameter r, we see that the 
spatial size of the traps increases with decreasing r with- 
out diverging. In the same way, the size of dynamical 



VI. CONCLUSION 

The coupled array of chaotic oscillators presented here 
has a number of properties that make it an important 
model. As was shown above and in Ref. 0, it is a faithful 
representation of the dynamics of an experimental set-up 
of coupled diode-resonators in chaotic regime character- 
ized by a stretched exponential distribution of trap time. 
Moreover, the low cost associated with solving numeri- 
cally the model of Eq. (||) allows us to study the system 
on time scale unreachable with atomic models; it takes 
less than a day, on a fast processor, to iterate a lattice of 
1000 sites over 10^ steps. With traps extending to more 
than 30,000 times the basic time step, long simulations 
are absolutely necessary to establish the nature of the 
dynamics in these systems. 

Results of these simulations demonstrate that the 
stretched exponential distributions arise from the compe- 
tition between a chaotic and a period-two regimes. The 
stretched-exponential requires some spatial organization 
to appear but does no imply diverging length scale as the 
traps become longer and longer: the period-4 structure 
is sufficient to stabilize a site onto a periodic orbit and 
the space of the spatial correlation is only weakly related 
to the length of the traps appearing in these systems. 

There are strong similarities between this system and 
the configurational glasses and we suggest that the un- 
derstanding gained here could be extended directly to 
these important materials. We are currently pursuing 
this avenue of research. 
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